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ABSTRACT 


A  theoretical  prediction  by  Dr.  Avihu  Ze’evi  of  the 
relative  contribution  to  the  optical  scintillation  by  different 
points  along  the  path  was  described  by  a  weight  function  for 
direct  and  exact  folded  path  spherical  wave  sources.  In  an 
effort  to  verify  this  prediction  a  turbulence  chamber  was  built 
to  allow  a  controlled  turbulence  source  to  be  moved  and 
measured  at  different  path  positions  in  conjunction  with 
scintillation  measurements.  The  experimental  results  follow 
Dr.  Ze'evi's  general  pattern  but  both  sources  are  less  weighted 
at  the  detector  end  than  predicted  and  the  folded  path  is  more 
heavily  weighted  at  the  target  end  than  predicted. 
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I .  INTRODUCTION 


A.  BACKGROUND 

In  optical  propagation  through  the  atmosphere  the  varying 
fluctuations  of  refractive  index  cause  random  deviations  in 
phase  across  a  propagation  wavefront.  This  change  of  the 
refractive  index  is  observed  as  the  "twinkling"  of  stars  or  the 
"image  boil"  of  lasers  and  is  known  as  scintillation. 
Scintillation  results  from  the  random  ray  bending  with  its 
associated  interference  manifesting  desultory  irradiance 
fluctuations  at  an  observation  plane. 

Theoretically,  in  the  weak  turbulence  region  an  accepted 
model  of  scintillation  for  spherical  waves  is  the  normalized 
log  amplitude  variance  CRef.  1]  given  by: 

a2  *  0.125  K7/,s  C2  L,l/6 
<  n 

where 

K  »  2ir/X ,  X  is  the  optical  wavelength 
L  *  Path  length  from  source  to  detector 

C2  *  Refractive  index  structure  constant  (a  direct 
n 

measure  of  the  level  of  atmospheric 
turbulence ) 

Experimentally,  scintillation  can  be  measured  by  placing  a 
point  detector  in  the  observation  plane  and  registering 
irradiance  level  fluctuations  with  time.  Scintillation 
strength  can  be  calculated  by  the  variance  of  the  normalized 
intensity  (irradiance),  as  follows: 
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<(I/Io)2>-  <I/Io>2 


I 

} 

°I/Io 

-  [<I2>-<I>2 ]/Io2 

-  Oj/Io 

where 

I  ■  I(r,t)  is  the  instantaneous  intensity  on  the  receiver 
a\  »  <I2>-<I>2 

IQ  »  <I(r,  t)>  is  the  ensemble  average 

Tatarski  CRef.  1]  found  that  in  regions  of  weak  turbulence  the 
log  irradiance  fluctuation  was  also  Gaussian  distributed.  Thus 
we  use  experimental  measurements  of  the  log  intensity  variance 
to  determine  the  scintillation  as  below* 

l  -  log  (I/l0) 

-  <12>-<1>2 

o[  •  <  (log  I/i0) 2>-<log  I/Iq>2 

*  <  (log  I* log  IQ)  2  >-<log  I  -  log  I Q>* 

*  <  (log  I  -  log  IQ)  (log  I  -  log  IQ)>  - 

<log  I  -  log  lQ>  <log  I  -  log  IQ> 

■  <  (log  I)2-  2  log  (IQ  +  I)  +  (log  IQ)  2>- 

[(<log  I>  ) 2-<2  log  IQ  log  I>+  (<log  IQ>)2] 

■  <(log  I}2>-  2  log  IQ  <log  I>  +  (log  IQ) 2  -<log  I>2 

+  2  log  Ic  <log  I>  -  (log  IQ) 2 
«  < (log  I)2>  -  <log  I>2 

a\  *  aLog  i 

This  allows  the  scintillation  to  be  described  by  the 
variance  of  log  intensity  fluctuations  received  at  a  detector 
with  no  dependence  on  the  average  signal  intensity  received. 
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Clifford  notes  [Ref.  2]  that  the  relationship  o£  =  4a*  is 
the  accepted  formula  to  correlate  the  empirical  (o^  )  and  the 
theoretical  ( aK  )  variances.  This  results  in  our  being  able  to 
describe  the  scintillation  for  spherical  waves  in  the  low 
turbulence  region  by: 

°‘l  =0.52K’/‘cn!  L"'* 

where  a2  ?  1  defines  the  low  turbulence  region. 

B.  PATH  WEIGHTING  THEORY 

The  works  of  Tatarski  [Ref.  1],  Fried  [Ref.  3  ],  Clifford 
[Ref.  2,  4] , Lutomirski  and  Yura  [Ref.  5],  Fante  [Ref.  6]  and 
others  have  led  to  a  fairly  clear  understanding  of  the 
processes  involved  for  scintillation  on  a  direct  path. 
Experimental  results  are  reasonably  close  to  theoretical 
predictions  for  direct  path  scintillation.  The  problem  of 
scintillation  on  folded  paths  is  less  understood,  and 
theoretical  predictions  of  existing  theories  have  not  been 
fully  confirmed  by  experimental  evidence. 

In  his  doctoral  dissertation  at  the  Naval  Postgraduate 
School,  Dr.  Avihu  Ze'evi  investigated  the  special  case  of 
optical  scintillation  on  folded  paths  displaced  a  small  angle  (9) 
from  the  direct  laser  path.  This  investigation  examined 
spherical  waves  and  a  plane  mirror  as  a  folding  target.  In 
particular,  Ze'evi  stated  that  for  the  direct  path  (one-way) 
laser  and  the  exact  (8*0)  folded  path  laser,  the  relative 


£ 


contribution  to  the  scintillation  of  different  points  along  the 
path  is  described  by  the  weight  function.  Ze'evi  stated  that 
the  weight  functions  W^(x)  for  spherical  waves  are: 

For  direct  path  WL(x)  =*  Cxd-x)]5^6 

For  the  folded  path  Wp(x)  =*  Cx(2-x)]5^6 

Figure  1  shows  the  optics  layout  for  the  direct  path  and  exact 
folded  path  lasers.  The  two  path -weighting  functions  are  shown 
in  Figure  2. 

*rom  Figure  2  we  can  see  that  in  the  case  of  the  exact 
folded  path  (8  =0),  the  main  contribution  to  the  scintillation 
comes  from  the  vicinity  of  the  mirror.  In  the  case  of  the 
single  path,  the  main  contribution  is  from  the  center  of  the 
path.  Ze'evi  further  noted  that  it  was  very  unlikely  that  we 
could  find  the  proper  conditions  under  which  C^(x)  along  an 
atmospheric  path  is  constant,  and  that  we  should  be  careful 
when  we  compare  results  from  two  types  of  path. 

The  development  and  use  of  a  controlled  environment  chamber 
to  test  Ze'evi's  predicted  scintillation  path  weighting 
functions  for  the  direct  path  and  exact  folded  path  spherical 
waves  was  the  main  goal  of  this  experiment. 
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II .  EXPERIMENTAL  APPROACH 

A.  GENERAL  CONSIDERATIONS 

As  a  result  of  the  continuing  experimental  studies  at  the 
Naval  Postgraduate  School  into  the  effect  of  turbulent 
atmosphere  on  optical  propagation  [Ref.  8]  and  the  attempts  by 
Hodgini  CRef.  9],  and  Speer  and  Parker  [Ref.  10]  to  validate 
portions  of  Ze'evi's  doctoral  dissertation  we  began  our 
analysis  with  access  to  experimental  apparatus  designed  to 
measure  scintillation  for  direct  path  and  folded  path  lasers. 

As  a  result  of  observing  the  efforts  of  Speer  and  Parker  we 
realized  a  controlled  environment  was  essential  to  produce 
meaningful  results  in  an  analysis  of  scintillation  path 
weighting.  The  turbulence  source  used  to  cause  scintillation 
along  our  path  needed  to  be  as  uniform  as  possible,  measurable 
for  C^.  in  order  to  normalize  our  results,  and  capable  of  being 
positioned  at  different  locations  along  the  laser  propagation 
path.  Additionally,  the  optics  needed  to  be  isolated  from  air 
currents  (wind)  that  had  caused  vibrations  in  past  experiments 
and  thus  introduced  added  uncertainties  to  the  results.  For 
these  reasons  we  decided  to  add  a  turbulence  chamber  to  the 
existing  apparatus.  This  chamber  consisted  of  a  "tunnel", 
a  controlled  turbulence  source  with  a  method  for  measuring  C-p 
of  the  disturbance,  and  protection  for  the  system's  optical 
components.  With  these  considerations  in  mind  we  will  now 
begin  an  analysis  of  the  entire  experimental  apparatus  to 
measure  the  scintillation  during  our  experiment. 
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B.  TURBULENCE  CHAMBER 


The  basic  concept  for  the  turbulence  chamber  called  for 
eliminating  turbulence  and  vibrations  to  optical  components 
unless  we  caused  the  turbulence  at  a  known  location  in  the 
system  and  could  measure  this  turbulence.  The  following 
components  comprised  our  turbulence  chamber: 

1.  Tunnel 

A  61  meter  long  tunnel  of  .61  meter  height  and  .61 
meter  width  was  constructed  by  joining  25  sections  of  2.44  m  by 
.61  m  by  .61  m  boxes  made  of  1/4  inch  plywood.  To  prepare  the 
tunnel  for  possible  outdoor  use  the  plywood  was  treated  with  a 
water  sealant  and  painted  white  on  the  exterior  to  lessen  the 
diurnal  heating  effects.  The  tunnel  interior  was  painted  flat 
black  to  eliminate  reflections.  Our  experiment,  however,  was 
conducted  in  the  145  meter  basement  of  Spanagel  Hall  thus 
rendering  some  of  these  precautions  unnecessary.  To  allow 
access  to  basement  laboratories  the  individual  boxes  were  only 
joined  together  to  form  a  tunnel  for  the  actual  experimental 
runs.  Box  seams  and  the  edges  that  mated  to  other  boxes  were 
sealed  with  76  mm  wide  insulation  tape  to  eliminate  unwanted 
turbulence  entrance.  One  tunnel  section  was  modified  to  house 
the  turbulence  source  and  will  be  discussed  later. 

2.  Turbulence  Source 

The  turbulence  source  was  chosen  by  Plenniken  [Ref.  11] 
to  insure  we  operated  in  an  isotropic  turbulence  region.  The 
turbulence  source  consisted  of  a  heat  element  connected  to  a 
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power stat  output  to  control  the  amount  of  heat  produced.  The 
heat  element  was  mounted  on  ceramic  insulators  below  a  screen 
with  .0015  meter  square  meshing  (Figure  3)  to  break  the 
turbulence  into  small  evenly  distributed  patterns.  At  a 
setting  of  70  volts  the  heater  produced  a  local  air  temperature 
of  27-30  Celsius  at  22  cm  above  the  screen. 

The  turbulence  source  was  mounted  in  a  modified  tunnel 
section  (Figure  4)  to  insure  its  configuration  would  remain 
essentially  constant  and  to  ease  in  relocating  the  turbulence 
source  during  experiments.  A  relatively  simple  procedure  was 
devised  where  this  section  could  be  entered  into  the  tunnel  at 
a  desired  location  and  the  tunnel  taped  together  thus  creating 
a  constant  tunnel  configuration  with  only  the  turbulence 
location  changing.  For  our  experiments  the  turbulence  source 
was  introduced  at  1.0,  5.8,  10.7,  15.5,  20.4,  25.3,  30.1,  35.0, 
39.8,  44.7,  49.5,  54.4,  and  60.0  meters  from  the  detector. 

This  change  of  turbulence  location  could  be  completed  in 
several  minutes  without  disturbing  or  adjusting  any  optical 
components  involved  in  the  experiment.  The  actual 
configuration  of  the  modified  tunnel  section  evolved  throughout 
this  experiment  as  a  result  of  our  desire  to  isolate  the 
turbulence  section  from  the  hallway  environmental  influences  as 
much  as  possible.  Modifications  will  be  discussed  in  the 
experimental  section. 

3.  Turbulence  Measurement  Equipment 

It  was  our  desire  to  have  the  turbulence  we  caused  to 
be  as  constant  as  feasible.  Realizing  the  random  nature  of  the 


Heat  Source  Showing  the  Heat 
Element  Mounted  Below  the  Screen 


Figure  4 


Optics  Protection  Cap  Joined  to  the 
Modified  Tunnel  Section  that  Housed 
the  Heating  Source.  The  platinum 
probes  for  Ci|  measurement  are  shown 
entering  the  tunnel  just  above  the 
baffle.  Heat  source  tunnel  section 
has  chimney  and  roof  mounted. 
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heat  dissipation  we  chose  to  measure  the  temperature  structure 
index  and  then  normalize  our  measured  experimental  variance 
for  the  fluctuations  of  C,|  at  the  different  path  locations. 

was  measured  by  inserting  2  platinum  probes  spared  6  cm 
apart  a  height  of  22  cm  above  the  heat  source  screen.  The 
platinum  probes  were  sampled  at  a  rate  of  11  Hertz  and  500 
samples  were  taken  to  compute  C^,  .  This  provided  a  45  second 
sampling  period  in  which  to  compute  an  accurate  picture  of  the 
temperature  structure  index.  Four  such  runs  were  made  and 
averaged  at  each  heat  source  location  just  prior  to  collecting 
scintillation  data.  Figure  5  shows  the  equipment  utilized  to 
process  the  data  from  the  platinum  probes.  For  details  on  the 
operation  the  equipment  utilized  to  calculate  see  Reference  11. 

4.  Optics  Protection 

To  attempt  to  eliminate  unwanted  turbulence  from  the 
tunnel  and  prevent  the  possibility  of  air  currets  vibrating  our 
optica  we  capped  the  tunnel  ends  (Figure  6).  Because  of  the 
size  of  our  optics  the  end  caps  had  to  be  larger  than  the  .6lm 
x  .61m  tunnel  section  so  we  chose  standard  laboratory  tables 
large  enough  to  house  the  optics  and  boxed  off  the  legs  with 
the  same  type  plywood  used  to  construct  the  tunnel.  These  were 
taped  in  an  effort  to  stop  all  turbulence  and  light  leaks. 

During  experimental  runs  these  end  caps  were  taped  to  the 
tunnel  and  all  data  and  power  wires  to  the  optics  were  run 
under  the  back  entrance  to  the  caps  and  tape  was  utilized  to 
seal  possible  turbulence  entrances. 
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Figure  5.  C^,  Calculation  Layout 
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Optics  Cap  for  Target  End  of  Tunnel.  Prom  left  direct  path 
laser,  polarizer  pair,  lens,  chopper  direct  path  beam  splitter, 
and  folded  path  target  mirror. 


Optics  Cap  for  Detector  End  of  Tunnel.  From  left,  detector, 
folded  path  beam  splitter,  polarizer,  Pockels  cell,  and 
laser.  Note  Pockels  modulation  system,  PAR,  and  high 
voltage  on  top  of  optics  cap. 

Figure  6.  Open  Optics  Caps 
21 


C.  OPTICAL  EQUIPMENT 


The  optical  setup  utilized  is  illustrated  in  Figure  1.  The 
target  end  consisted  of  the  direct  path  laser,  whose  intensity 
was  controlled  by  use  of  two  polarizers.  Next,  the  direct  path 
beam  was  diverged  by  a  lens,  and  modulatd  by  a  mechanical 
chopper.  The  modulated  beam  was  then  reflected  down  the  tunnel 
to  the  detector  aperture  by  a  beam  splitter.  Located  directly 
behind  the  beam  splitter  was  the  target  mirror  for  the  folded 
path.  The  detector  end  optics  consisted  of  the  folded  path 
laser,  Pockels  cell  modulator,  a  polarizer,  a  beam  splitter  to 
cause  the  direct  path  and  folded  path  beams  to  be  coincident  as 
the  folded  path  beam  travelled  the  tunnel,  passed  through  the 
direct  path  splitter,  reflected  off  the  target  mirror  and 
returned  down  the  tunnel  on  the  same  path  to  the  detector  as 
the  direct  path  laser.  The  aperture  and  detector  were  located 
behind  the  folded  path  splitter. 

1 .  Laser  Sources 

The  two  lasers  used  throughout  the  experiment  were 
HeNe  lasers  with  a  wavelength  of  .6328  micrometers.  The  low 
output  of  these  lasers,  .95  milliwatt,  made  them  eyesafe  and 
thus  reduced  the  safety  hazard  and  associated  problems  in 
both  alignment  and  operation.  Both  lasers  were  mounted  on 
platforms  controlled  by  micrometer  screw  adjustments  to 
facilitate  precision  alignment.  It  was  essential  to  utilize 
laser  sources  whose  internal  stability  would  be  such  as  not  to 
effect  the  scintillation  we  were  measuring.  Hodgini  measured 


the  source  variances  of  these  lasers  as  less  than  3  x  10 


and  ratio  of  the  source  variance  of  these  two  lasers  was  .996 
[Ref.  9].  Although  the  lasers  were  noted  to  become  unstable 
with  time*  the  correlation  of  the  two  sources  when  chopped 
alternately  with  a  pulse  separation  of  0.3  milliseconds  over 
the  same  optical  path  did  not,  and  thus  they  could  be  used 
together  in  measuring  scintillation  with  some  confidence  that 
their  internal  instability  would  not  contribute  to  the 
scintillation  we  were  measuring. 

2.  Target  End  Optics 

The  target  end  optics  consisted  of  a  set  of  adjustable 
polarizers  for  fine  signal  intensity  reductions  to  match  the 
direct  path  intensity  to  that  of  the  folded  path  laser.  A  lens 
was  utilized  to  expand  the  beam  to  prevent  beam  wander  from 
affecting  the  scintillation  measurements.  The  four-inch  beam 
splitter  was  a  half-inch  thick  optical  window  with  an 
advertised  flatness  correct  to  1/10  of  a  wavelength  and 
parallelism  within  two  seconds  of  arc.  An  antireflective 
coating  was  applied  and  one  side  of  the  splitter  silvered  to 
allow  50  percent  transmission.  The  target  mirror  for  the 
folded  path  laser  was  a  five  inch  optical  flat  of  good  optical 
flatness. 

3.  Detector  End  Optics 

The  folded  path  beam  passed  immediately  through  a 
coherent  Associates  Model  3003  Modulation  System  Pockels  cell. 
This  Pockels  cell  in  conjunction  with  the  polarizer  that  the 
beam  passed  through  next  were  used  to  modulate  the  laser  beam. 


D.  PULSE  FORMING  AND  DETECTION  EQUIPMENT 

The  pulse  forming  equipment  in  our  experiment  included  a 
mechanical  chopper  with  an  open  to  close  ratio  of  1:4  that 
rotated  at  a  3000  rpm  rate.  This  provided  a  chopping  frequency 
of  1000  Hertz.  The  chopper  rotation  modulated  a  light  emitting 
diode  mounted  opposite  a  detector.  The  detector  sent  a  trigger 
signal  sent  to  a  Pulse/Function  Generator  thus  establishing 
the  timing  for  the  entire  system.  This  trigger  signal  was  then 
delayed  for  .4  milliseconds  before  triggering  the  Pockels  Cell 
Modulation  System  for  the  folded  path  laser.  This  timing 
allowed  the  direct  path  and  folded  path  lasers  to  arrive 
sequentially  at  the  detector  aperture.  The  aperture  size  was  3 
millimeters  and  this  was  well  within  the  6.2  millimeters 
required  to  be  within  the  first  Fresnel  Zone.  Zone  size  was 
calculated  as  follows: 
lF  =  (al)1/2 

£p  is  the  Fresnel  zone  radius 
X  =  .6328  micrometers 
L  =  61  meters 

By  utlizing  an  aperture  size  smaller  than  the  Fresnel  Zone 
and  smaller  than  the  10.2  mm  lateral  coherency  length  for  our 
HeNe  lasers  we  approximated  our  detector  as  a  point  detector  in 
the  observation  plane.  For  a  more  detailed  discussion  of  the 
choice  of  a  3  millimeter  aperture  see  Reference  9. 

Signals  passing  through  the  aperture  were  focused  on  the 
detector.  The  detector  was  an  RCA  silicon  photodiode  type 


C30872.  When  operated  at  370  volts  this  detector  provided  an 
extremely  high  signal  to  noise  ratio  without  the  need  for 
cooling. 

Next  the  output  of  the  detector  was  amplified  by  a 
Princeton  Applied  Research  (PAR)  Model  113  Low  Noise  Amplifier 
capable  of  yielding  a  signal  gain  of  up  to  10000.  Additionally, 
the  PAR  had  an  adjustable  bandwidth  with  high  and  low  frequency 
roll-offs.  The  roll-offs  were  carefully  set  to  eliminate  system 
noise  without  distorting  the  signal  sent  to  the  processing 
equipment.  Figure  7  gives  a  schematic  layout  of  the  pulse 
forming  and  detection  equipment  plus  the  respective  signal 
channeling  to  the  signal  processing  and  data  reduction 
equipment . 

E.  SIGNAL  PROCESSING  EQUIPMENT 

The  PAR  output  signal  was  sent  directly  to  the  direct  path 
and  folded  paths  demodulators.  A  trigger  pulse  activated  the 
demodulator  to  demodulate  only  its  respective  input  signal. 

The  demodulator  utilized  a  sample  and  hold  to  take  the 
instantaneous  maximum  of  the  incoming  signal.  The  background 
signal  was  subtracted  from  the  instantaneous  maximum  and  this 
difference  was  output  from  the  demodulator.  This  differential 
output  meant  that  the  true  maximum  signal  intensity  was 
available  for  further  processing.  Figure  8  shows  the  signal 
forms  and  timing  sequence  of  the  demodulator  waveforms. 
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Figure  7.  Diagram  of  Pulse  Forming,  Processing 
and  Data  Reduction  Equipment 


The  differential  output  of  the  demodulators  was  sent  to  the 
direct  and  folded  path's  respective  Hewlett-Packard  Model  7562A 
Log  Converters  since  we  desired  to  have  the  log  intensity 
variance  as  the  actual  measure  of  the  scintillation. 

Finally,  these  processed  signals  had  to  be  fed  to  the  data 
reduction  system  in  order  to  calculate  the  scintillation.  This 
process  was  complicated  by  the  fact  that  the  Nicolet  1080 
(NIC80)  Stored  Program  Computer  had  only  one  analog  to  digital 
converter  port.  To  time  share  this  port  for  use  by  both  laser 
signals  the  Atmospheric  Physics  Group  at  the  Naval  Postgraduate 
School  had  built  an  electronic  switch  which  allowed  the  Nicolet 
1080  to  accept  both  the  direct  path  and  folded  path  signals. 
Details  of  this  switch  operator  are  available  in  Reference  10. 
Figure  9  shows  the  Signal  Processing  Equipment  layout. 

F.  DATA  REDUCTION  EQUIPMENT 

The  computer  system  used  to  analyze  the  signal  was  composed 
of  four  components:  The  Nicolet  Instrument  Corporation  Model 

NIC-80;  the  Hewlett  Packard  (HP)  9825B  Caclulator;  the  HP-9871A 
Printer;  and  the  HP9862A  Calculator  Plotter.  Figure  10  shows 
the  Data  Reduction  Equipment  setup. 

The  program  for  the  NIC-80  is  available  in  Reference  10. 
During  experimental  runs  the  NIC-80  built  a  histogram  of  the 
number  of  counts  versus  point  signal  intensity  for  each  laser. 
This  was  accomplished  for  four  separate  runs  at  each  heater 
location.  Each  run  consisted  of  16,384  samples  per  laser 
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Figure  9.  Signal  Processing  Equipment 


Left  bay,  top  to  bottom:  Direct  path  demodulator, 

direct  path  log  converter, 
and  wavetek  signal  generator. 

Right  bay,  top  to  bottom:  Electronic  switch,  folded 

path  demodulator,  folded  path 
log  converter,  P-12's, 
oscilliscope,  and  P-25. 


Figure  10.  Data  Reduction  Equipment 


Front  row  from  left:  HP  987 1A  printer,  HP  98253 

calculator,  and  HP  9862A 
calculator  plotter. 
Background,  Nicolet  1080. 
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source  .  At  the  end  of  the  four  rune  this  data  was  transferred 
to  magnetic  tapes  loaded  in  the  HP-982 5B  for  storage  and 
statistical  analysis  upon  operator  command. 

The  HP-982 5B  calculator  was  the  system  controller  for  the 
NIC-80,  the  printer,  and  the  plotter.  Control  of  the  NIC-80 
was  achieved  through  a  computer  interface  designed  by  Professor 
Edmund  Milne  at  NPS.  The  HP9825-B  scintillation  program  is 
presented  in  Reference  10. 

To  calculate  the  statistical  parameters  we  desired  the 
intensity  distributions  from  each  run  were  loaded  from  their 
magnetic  tape  storage  locations  and  the  HP-9825B  processed  them 
in  accordance  with  accepted  formulas  for  mean,  variance, 
standard  deviation,  and  chi-square.  Speer  [Ref.  10]  determined 
the  for  a  95%  confidence  interval  the  chi-square  calculation 
with  a  goodness  of  fit  parameter  of  288.87  or  less  was  necessary 
to  assume  our  log  normal  curve  corresponded  to  a  true  Gaussian 
curve. 

Upon  completion  of  the  statistics  calculation  the  HP98253 
automatically  caused  the  printer  to  list  header  information 
detailing  the  date  the  experiment  was  conducted,  the  fact  that 
the  aperture  was  in  place  (AP),  the  detector  avalanche  voltage, 
PAR  bandpass,  and  the  fact  the  tests  were  for  the  exact  folded 
path  (0).  All  runs  were  conducted  with  a  PAR  gain  of  500. 
Additionally,  the  time  the  runs  began  were  noted  as  well  as  the 
turbulence  source  distance  from  the  detector.  In  the  columns 
below  this  are  listed  the  number  of  data  points  that  actually 


comprise  the  historgram,  the  variance,  and  the  chi-square 
value.  Four  runs  were  conducted  at  each  location  and  in  each 
case  the  direct  path  laser  parameters  are  listed  first,  then 
the  folded  path  parameters. 

Upon  completion  of  the  statistics  printout  of  each  run,  the 
HP9825B  presented  the  operator  a  choice  of  an  immediate  plot  of 
the  data  or  a  continuation  of  the  statistics  calculations .  The 
plotter  graphically  displayed  the  actual  point  by  point 
histogram  that  the  NIC-80  had  built  then  drew  a  solid  line  to 
show  a  true  normal  curve  calculated  by  the  chi-square  approach. 
Samples  of  the  statistical  printouts  for  these  experiments  are 
presented  in  Appendix  A  through  D.  These  sample  printouts  are 
accompanied  by  a  typical  direct  path  and  folded  path  plot  for 
that  turbulence  position.  Again,  each  plot  includes  a  header 
code  as  described  for  the  printout. 
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III.  EXPERIMENTAL  WORK 


A.  GOALS 

There  were  two  main  goals  for  this  experiment.  First,  we 
wanted  to  add  a  component  to  the  existing  scintillation 
measuring  system  that  would  allow  the  introduction  of  known 
turbulences  along  the  propagation  path  in  order  that  the 
Atmospheric  Physics  Research  Group  at  the  Na^al  Postgraduate 
School  could  investigate  various  optical  propagation  phenomena 
requiring  such  a  component.  Secondly,  we  wanted  to  compare  our 
experimental  results  for  direct  path  and  folded  path 
scintillation  weighting  with  the  theoretical  predictions  of 
Ze'evi. 

B .  PREPARATION 

In  order  to  insure  proper  equipment  procedures  and  to 
exercise  experimental  techniques,  numerous  runs  were  conducted 
utilizing  the  overhead  145  m  range  in  Spanagel  Hall.  This 
range  is  below  eight  overhead  ventillation  shafts  and  is 
adjoined  by  many  hallway  doorways  and  four  building  entrances. 
Thus,  there  is  a  great  deal  of  turbulence  present  but  the 
turbulence  was  not  controllable.  No  data  from  these  runs  are 
included  in  the  report. 

Next  we  transferred  the  optics  portions  of  the  experiment 
to  the  basement  floor  location  of  Spanagel  where  the  61  meter 
tunnel  would  be  utilized.  Practice  runs  were  conducted  during 
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daylight  working  hours  with  the  tunnel  sections  open.  Again 
the  turbulence  was  uncontrollable  but  our  runs  showed  direct 
path  variances  in  the  .014  to  .068  range  and  folded  path 
variances  in  the  range  of  .018  to  .118.  These  are  mentioned 
only  to  give  some  basis  by  which  to  compare  the  quieting  effect 
of  our  tunnel  discussed  next.  Prior  to  each  experiment  we 
sealed  the  entire  tunnel,  sealed  all  adjoining  doorways,  and 
measured  a  “quiet"  tunnel  reading  that  we  utilized  as  our 
background  reading  for  that  experiment.  The  average  quiet 
tunnel  variances  for  all  four  experiments  were  .003  for  the 
direct  path  laser  and  .006  for  the  folded  path  laser  thus 
giving  us  reasonable  confidence  that  our  tunnel  was  indeed  a 
region  that  controlled  the  turbulence  the  lasier  had  to 
propagate  through.  By  subtracting  each  experiment ’ s  respective 
“quiet"  tunnel  variance  from  the  individual  turbulence  position 
variance  during  actual  runs,  we  sought  to  achieve  the  true 
effect  of  the  turbulence  at  that  particular  position.  With 
these  considerations  in  mind  let  us  now  discuss  our  four 
experimental  runs. 

C.  FIRST  EXPERIMENT  -  16  MARCH  1983 

For  the  16  March  1983  experiment,  the  equipment  setup  as 
described  in  paragraph  II  with  the  following  noteworthy  ponts . 
The  lens  utilized  in  the  direct  path  setup  was  a  285  mm  plane 
convex  lens  set  a  focal  length  from  the  chopper  wheel.  The 
turbulence  source  tunnel  section  was  missing  the  chimney,  roof, 
and  baffle  shown  in  Figure  4.  Throughout  the  run  we 
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attempted  to  vary  the  voltage  applied  to  the  heating  coil  to 
stabilize  .  Voltages  utilized  varied  from  90  to  105  volts 
but  we  failed  to  truly  stabilize  C^,  .  As  a  result  of  this 
problem,  we  decided  to  normalize  all  data  for  Cij,  changes. 
Figures  11  through  14  show  the  positional  variances  versus 
turbulence  distance  from  the  detector  for  this  experiment.  The 
solid  lines  on  each  graph  are  the  least  squares  fit  of  Ze'evi's 
theoretical  path  weighting  curves.  These  plot  parameters  apply 
to  the  plots  accompanying  all  four  experiments.  Detailed  data 
for  this  run  is  located  in  Appendix  A. 

D.  SECOND  EXPERIMENT  -  27  MARCH  1983 

In  analyzing  the  data  from  the  16  March  1983  run  we  decided 
to  add  the  chimney  to  the  tunnel  to  further  isolate  the  tunnel 
and  turbulence  source  from  the  possibility  that  hallway  breezes 
may  have  introduced  fluctuations  in  the  turbulence  presented  to 
the  lasers  during  propagation.  Additionally,  we  realized  the 
interference  patterns  we  observed  on  the  detector  face  for  the 
direct  path  and  folded  path  patterns  may  have  influenced  our 
readings.  To  eliminate  these  interference  patterns,  great  care 
was  taken  to  eliminate  the  reflections  off  the  back  face  of  the 
splitters.  The  folded  path  interference  pattern  was  eliminated 
by  intercepting  the  secondary  beam  off  the  back  face  of  the 
folded  path  splitter  with  a  small  wooden  stick.  Careful 
examination  of  the  beam  structure  at  the  target  mirror  and  on 
the  detector  face  indicated  that  this  eliminated  the 
interference  problem  without  introducing  structure  to  the  beam. 


Figure  14.  Plot  of  the  Normalized  Folded  Path  Variance  Versus  Disturbance 
Distance  from  the  Detector  for  the  16  March  1983  Experiment. 
Data  was  normalized  for  Ct*  Solid  line  represents  the  least 
square  fit  of  Ze’evi's  prediction. 


The  direct  path  pattern  was  eliminated  by  removing  the  285  mm 
lens  and  properly  orienting  the  polarizers  and  direct  oath 
splitter  to  achieve  the  Brewster  angle.  This  eliminated  the 
interference  pattern  but  in  analyzing  the  results  of  this 
experiment  we  realized  that  we  had  introduced  and  measured  beam 
wander  as  a  result  of  too  narrow  a  beam  in  the  portion  of  the 
tunnel  near  the  direct  path  splitter.  However,  the  folded  path 
data  for  the  run  was  valid  and  the  variance  and  normalized 
variance  plots  are  presented  in  Figures  15  and  16.  Detailed 
data  for  the  27  March  1983  run  is  presented  in  Appendix  B. 

E.  THIRD  EXPERIMENT  -  30  MARCH  1983 

To  eliminate  the  beam  wander  and  interference  pattern 
problem  for  the  direct  path  laser  we  introduced  a  10  cm  focal 
length  convex-convex  lens  as  our  lens  before  the  chopper.  This 
gave  us  a  large  uniformly  structured  beam  at  the  aperture  for 
our  third  experiment.  Additionally  we  decided  to  lower  the 
heater  voltage  to  70  volts  for  this  run  and  to  leave  the 
voltage  constant  throughout  the  experiment.  The  plots  of  the 
variances  versus  distrubance  distance  for  this  experiment  are 
presented  in  Figures  17  through  20.  Detailed  data  for  this  run 
is  located  in  Appendix  C. 

F.  FOURTH  EXPERIMENT  -  13  APRIL  1983 

In  analyzing  the  data  of  30  March  1983  we  decided  to  add 
the  roof  and  baffle  to  the  turbulence  source  tunnel  to  better 
isolate  the  turbulence  source  from  the  hallway  environment. 
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Figure  16.  Plot  of  Normalized  Folded  Path  (two-way)  Variance  Versus  Disturbance 
Distance  from  the  Detector  for  the  27  March  1983  Experiment. 

Data  was  normalized  for  C^.  Solid  line  represents  the  least 
squares  fit  of  Ze'evi’s  prediction. 


jure  18.  Plot  of  Normalized  Direct  Path  (one-way)  Variance  Versus  Dis' 
Distance  from  the  Detector  for  the  30  March  1983  Experiment. 
Data  was  normalized  for  C,p.  Solid  line  represents  the  least 
fit  of  Ze'evi's  prediction. 
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We  feared  that  air  from  an  overhead  vent  or  cooler  air  from 


building  stairwells  may  have  directly  entered  the  heater. 
Additionally,  realizing  that  the  results  we  had  received  thus 
far  showed  a  consistent  pattern  of  less  weighting  at  the 
detector  end  for  both  sources  and  greater  weighting  at  the 
target  end  for  the  folded  path  laser  than  predicted  by  Ze’evi 
we  decided  to  reverse  the  optics  locations  in  the  hallway  to 
insure  that  our  results  were  not  dependent  on  a  particular 
hallway  configuration.  The  plots  of  variances  versus 
disturbance  location  from  the  detector  are  presented  in  Figures 
21  through  24.  It  is  significant  to  note  that  reversing  the 
optics  did  not  change  the  overall  path  weighting  pattern  we 
previously  encountered.  Detailed  statistics  of  the  13  April 
1983  run  are  located  in  Appendix  D. 
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Figure  22.  Plot  of  Normalized  Direct  Path  (one-way)  Variance  Versus  Disturbance 
Distance  from  the  Detector  for  the  13  April  1983  Experiment. 

Data  was  normalized  for  C<f>.  Solid  line  represents  the  least 
squares  fit  of  Ze'evi's  prediction. 


Plot  of  Folded  Path  (two-way)  Variance  Versus  Disturbance 
Distance  from  the  Detector  for  the  13  April  1983  Experiment. 

Solid  line  represents  the  least  squares  fit  of  Ze'evi's  prediction 


IV.  CONCLUSIONS 


A.  FINDINGS 

Based  upon  our  experimental  results  and  data  analysis  we 
state  the  following  findings: 

1.  In  examining  our  four  experimental  runs  we  conclude  that 
the  data  we  obtained  follows  the  general  pattern  of  Ze’evi's 
prediction  but  fails  to  conform  exactly  to  his  predicted  path 
weighting.  In  particular,  our  data  indicates  a  lower  path 
weighting  in  the  vicinity  of  the  detector  for  both  the  direct 
path  and  folded  path  and  that  a  much  more  severe  path  weighting 
in  the  vicinity  of  the  target  mirror  existed  than  Ze'evi 
predicted. 

2.  Our  experiments  were  conducted  in  the  low  turbulence 
region  defined  by  <  1  .  The  actual  o ^  measured  were 
generally  less  than  .9.  In  examining  our  chi-square  values  for 
these  distributions  we  must  conclude  that  in  the  turbulence 
region  of  our  experiment  the  log  intensity  distribution  was  not 
a  true  Gaussian  distribution  since  our  chi-square  goodness  of 
fit  parameters  generally  exceeded  the  established  value  of 
288.87  for  a  95%  confidence  interval. 

3.  We  state  with  confidence  that  the  turbulence  chamber 
setup  utilized  was  an  acceptable  method  for  introducing  a  known 
turbulence  in  the  laser  propagation  paths  and  that  run  number 
four  demonstrated  that  our  findings  are  independent  of  the 
terminal  optics  locations  with  respect  to  the  hallway 
configuration. 
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4.  A  closer  examination  of  the  folded  path  variance  versus 
disturbance  distance  indicated  that  a  complicated  three  region 
pattern  may  actually  exist.  Thin  may  be  worthy  of  further 
examination. 

5.  As  we  attempted  to  refine  our  experiment  and  believed  we 
understood  more  about  what  was  actually  happening  along  the 
propagation  path,  we  ^ecran  to  fully  realize  that  unpredictable 
scintillation  effects  could  occur  at  any  time.  As  an  example 
we  could  not  explain  the  unusual  dip  in  the  direct  path 
variance  for  the  30.1  meter  position  in  our  last  run. 

B.  FUTURE  WORK 

Future  work  utilizing  the  system  as  currently  configured 
could  include: 

1.  Repetition  of  the  current  experiment  to  build  a  large 
data  base  may  be  considered  in  order  to  aid  in  ascertaining  a 
path  weightina  formula  that  more  accurately  predicts  the 
empirical  data. 

2.  The  use  of  the  turbulence  chamber  to  attempt  to  verify 
Ze’evi's  predicted  correlation  of  direct  and  folded  path 
scintillation  for  small  ancles  other  than  the  exact  folded  path 
(9=0)  would  be  of  interest. 

3.  Use  of  a  corner  cube  reflector  and  other  specular 
reflectors  for  the  folded  path  target  and  the  examination  of 
the  scintillation  path  weightina  for  such  setups  would  be  of 
value  to  many  currently  applied  uses  of  lasers  utilizing  folded 
paths  . 
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Figure  25.  Scintillation  Statistics 
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Figure  26.  Scintillation  Data  and  Theoretical  Curve 
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Scintillation  Data  and  Theoretical  Curve 


Scintillation  Data  and  Theoretical  Curve 
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Figure  32.  Scintillation  Data  and  Theoretical  Curve 
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Figure  33.  Scintillation  Data  and  Theoretical  Curve 
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Figure  38.  Scintillation  Data  and  Theoretical  Curve 


Figure  39.  Scintillation  Data  and  Theoretical  Curve 
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Figure  40.  Scintillation  Statistics 


Scintillation  Data  and  Theoretical  Curve 


Figure  42.  Scintillation  Data  and  Theoretical  Curve 
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Figure  43.  Scintillation  Statistics 


Figure  44.  Scintillation  Data  and  Theoretical  Curve 
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Scintillation  Data  and  Theoretical  Curve 
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Figure  50.  Scintillation  Data  and  Theoretical  Curve 


nor  or  2  w»r 
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Figure  54.  Scintillation  Data  and  Theoretical  Curve 
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